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DIFFUSIVE LIMIT WITH GEOMETRIC CORRECTION OF UNSTEADY NEUTRON 

TRANSPORT EQUATION 


LEI WU 


Abstract. We consider the diffusive limit of an unsteady neutron transport equation in a two-dimensional 
plate with one-speed velocity. We show the solution can be approximated by the sum of interior solution, 
initial layer, and boundary layer with geometric correction. Also, we construct a counterexample to the 
classical theory in [I] which states the behavior of solution near boundary can be described by the Knudsen 
layer derived from the Milne problem. 

Keywords: compatibility condition, e-Milne problem, Knudsen layer, geometric correction. 


1. Introduction and Notation 


1.1. Problem Formulation. We consider a homogeneous isotropic unsteady neutron transport equation 
in a two-dimensional unit plate 11 = {af = (xi,X 2 ) '■ |af| < 1} with one-speed velocity E = {w = (wi,W 2 ) '■ 
w G 5^} as 



0 in [0, oo) X n, 
h{x, w) in n 

g(t,xo,w) for w ■ n < 0 and xq G dfl, 


( 1 . 1 ) 


where 



( 1 . 2 ) 


and n is the outward normal vector on dfl, with the Knudsen number 0 < e << 1. The initial and boundary 
data satisfy the compatibility condition 


(1.3) 


h{xo,w) = g{0, xo,w) for w-n<0 and xq € dfl. 


We intend to study the diffusive limit of as e —?> 0. 


Based on the flow direction, we can divide the boundary F = {(a?, ■ul) : x G 911} into the in-flow boundary 
F”, the out-flow boundary F+, and the grazing set F° as 

(1.4) F“ = {(af, rZl): x G dV,, w ■ n < 0}, 

(1.5) F’*' = {(ai, fU) : X G 911, w ■ n> 0}, 

(1.6) F° = {(x, iZl): X G 911, w; • fi = 0}. 

It is easy to see F = F+ U F“ U F°. 


The study of neutron transport equation dates back to 1950s. The main methods include the explicit 


formula and spectral analysis of the transport operators(see [5], [4], [^, [7], [8], [9], [10], [TT], [H])- In the 
classical paper [1], a systematic construction of boundary layer was provided via Milne problem. However, 
this construction was proved to be problematic for steady equation in [13] and a new boundary layer con¬ 
struction based on e-Milne problem with geometric correction was presented. In this paper, we extend this 
result to unsteady equation and consider a more complicated case with initial layer involved. 

1.2. Main Results. We first present the well-posedness of the equation (HU). 

Theorem 1.1. Assume g{t,Xo,w) G L°°([0, oo) x F“) and h{x,w) G x S^). Then for the unsteady 

neutron transport equation im, there exists a unique solution u^{t,x,w) G L°“([0,oo) x fl x S^) satisfying 


(1.7) 
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Then we can show the diffusive limit of the equation (HH). 

Theorem 1.2. Assume g{t,xo,w) G C'"‘([0,c») x F”) and h{x,w) G x 5^). Then for the unsteady 

neutron transport equation C3D, the unique solution u’^(t,x,w) G L°°([0,(X)) x x 5^) satisfies 

(1-8) ||u" - t/o - - '^B,o|li=o = o(l)> 

where the interior solution Uq is defined in US. 56\) . the initial layer is defined in i3.55\) . and the boundary 
layer g is defined in \3.54^ . Moreover, if g(t, 9, (f) = f^e~^ cos (f and h(x, w) = 0, then there exists aC > 0 
such that 

(1.9) llu*^ — Uq — — '^b,oILoo > C > 0 

when e is sufficiently small, where the interior solution Uq is defined in i6.19\) . the initial layer is defined 
in h6.18\) . and the boundary layer f/g q is defined in Jg. 17] ). 

Remark 1.3. 9 and 4> are defined in and VS. 39\) . 

It is easy to see, by a similar argument, the results in Theorem 11.11 and Theorem 11.21 also hold for the 
one-dimensional unsteady neutron transport equation, where the temporal domain is [0,oo), spacial domain 
is [0, L] for fixed L > 0, and velocity domain is [—1/2,1/2]. 

1.3. Notation and Structure of This Paper. Throughout this paper, C > 0 denotes a constant that 
only depends on the parameter Q, but does not depend on the data. It is referred as universal and can 
change from one inequality to another. When we write C{z), it means a certain positive constant depending 
on the quantity 0 . We write a < b to denote a < Cb. 

Our paper is organized as follows: in Section 2, we establish the L°° well-posedness of the equation CID 
and prove Theorem ll.il in Section 3, we present the asymptotic analysis of the equation (HI; in Section 4, 
we give the main results of the e-Milne problem with geometric correction; in Section 5, we prove the first 
part of Theorem 1 1.21 finally, in Section 6, we prove the second part of Theorem 11.21 


( 2 . 1 ) 


2. Well-posedness of Unsteady Neutron Transport Equation 
In this section, we consider the well-posedness of the unsteady neutron transport equation 
e^dtu + ew ■'VxU + u — u = f{t,x,w) in [0,oo)xr2, 
u(0,x,uJ) = h(x,w) in U 

u(t, xo,w) = g(t,xo,w) for w ■ n < 0 and xq G dfl, 
where the initial and boundary data satisfy the compatibility condition 

(2.2) h{xo,w) = g{0,Xo,w) for w-h<0 and xq G dV,. 

We define the and L°° norms in U x 5^ as usual: 

\ 1/2 

(2.3) ll/llL2(nx5i) 

(2.4) 


[ [ |/(f,w)|^du;df) , 

Jn Js^ J 


lL”(nx5i) 


sup \f{x,w)\. 
{x,w)^QxS^ 


Define the and L°° norms on the boundary as follows: 


(2.5) 

1/ lL2(r) 

(2.6) 

1/ lL2(r±) 

(2.7) 

1/ lL“(r) 

(2.8) 

1/ lBoo(r±) 


JJ |/(a:, ru)]^ luJ • n| duJdT^ , 

JJ \f{x,w)f\w-h\dwdx 


1/2 

1/2 


= sup |/(x,u;)|, 


Similar notation also applies to the space [0,oo) x U x 5^, [0,oo) x F, and [0, oo) x F=*=. 
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2.1. Preliminaries. In order to show the and L°° well-posedness of the equation (EH), we start with 
some preparations of the penalized neutron transport equation. 

Lemma 2.1. Assume f{t, x, w) € oo) x 11 x 5^), h{x, w) G L°°(f2 x 5^) and g{t, Xq, w) G L°°([0, oo) x 

r“). Then for the penalized transport equation 

{ Xux + e'^dtUx + ew -VxUx + u\ = f{t,x,w) in [0, oo) x H, 
ux{0,x,w) = h{x,w) in 12 

ux(t,Xo,w) = g{t,Xo,w) for w ■ n < 0 and xq G dV,. 
with A > 0 as a penalty parameter, there exists a solution ux{t,x,w) G L°°([0,T] x 11 x 5^) satisfying 
(2.10) l|■“A|lL~([0,oo)xnxSl) — ll5llLoo([o,oo)xr-) + ll^llL“(nxSi) + ll/llL“([0,oo)xnx5i) • 

Proof. The characteristics (T(s), X(s), lT(s)) of the equation (12.91) which goes through (f, x, w) is defined by 

(T(0),X(0),1T(0)) = it,x,w) 

dT{s) 2 


( 2 . 11 ) 


which implies 


ds 

dX{s) 

ds 

dW{s) 

ds 


= elT(s), 
= 0 . 



r T{s) 

= t + e^s, 

(2.12) 

{ X{s) 

= X + (ew)s 


[ Vk(s) 

= w, 


Hence, we can rewrite the equation EH) along the characteristics as 
(2.13) ux{t,x,w) 


= 1 


{t>e 2 tt}(^ 9 (t - - etbw,w)e + J f(^t-e‘^{tb-s),x-e{tb-s)w,w)e *)ds^ 

+I{t<e 2 tfc} (etw)/e^,'(fi)e“^^+^^‘/'' + J /(e^s, f - e(t/e^ - s)rZ;, , 


where the backward exit time tb is defined as 

(2.14) tb{x,w) = inf{s > 0 : (a; — esw,w) G T”}. 

Then we can naturally estimate 

(2-15) II^A|lioo([Q oo)xnxsi) 






Il5ll 


L“([0,oo)xr- 


+I{t<62t6}^e A+^)t/^ ll^llL=(nx5i) - 1 


^ 111 llL°°([0,oo)xnx5i) 

\ _ 0 ( 1 +-^)*/ 

ll/llL°°([0,oo)xnx5i) 


^ ll5'llL°°([o,oo)xr-) + ll^llL“(nx5i) + ll/llL“([o,oo)xnxSi) • 

Since ux can be explicitly traced back to the initial or boundary data, the existence naturally follows from 
above estimate. □ 

Lemma 2.2. Assume f{t, x, w) G L°“([0, 00 ) x fl x 5^), h{x, w) G x S^) and g{t, Xq, w) G L°“([0, 00 ) x 

T”). Then for the penalized neutron transport equation 

Xux + e^dtux + ew ■'S/xUx + ux — ux = f{t,x,w) in [0, oo)xll, 

(2.16) ux{0,x,w) = h{x,w) in H 

ux{t,xo,w) = g{t,xo,w) for xq G 911 andw-n<0. 
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with A > 0^ there exists a solution ux{t,x^w) G L°°([0,cx)) x ^ x S^) satisfying 

(2-17) hAll L“([0,oo)xnx5i) — ll5llLoo([0,oo)xr-) + ll^llL“(nx5i) + ll/llL“([0,oo)xnx5i) ^ • 

Proof. We define an approximating sequence {u^}'^q, where = 0 and 

Xu^l + e'^dtu’l + ew ■ Vxu'l + u'l — = f(t,x,w) in [0,oo)xri, 

u'^{0,x,w) = h{x,w) in V. 


(2.18) 


l(t,xo,w) = g(t,xo,w) for xq € dfl andw-n<0. 


By Lemma [2.11 this sequence is well-defined and llu* 


A|lL=°([0,oo)xnxSi) 


< oo. 


The characteristics and the backward exit time are defined as (12.111) and (I2.14L so we rewrite equation 
(12.181) along the characteristics as 

(2.19) 

u\(t,x,w) 

= -b J + f){t - e^(tb - s),x - e(tb - 

+l[t<eH,}(h(x- (etu()/e2,u;)e-(i+^)‘/^^ x - e(t/e^ - 

We define the difference = u\ — for /c > 1. Then satisfies 

(2.20) v’""^'^(x,w) = l{t>eHb}(^J Vx~^it - e‘^(tb - s),x - e(tb - 


+ I{t<e2t^} 


t/e 


;r,k-lr^2 


(eh, X - e(t/e^ - s)w, ""Ms 


Since II^ooqq oo)xnxSi) - Ih IIl“([o,oo)x 


( 2 . 21 ) 


0 x 51 )’ directly estimate 

^max{t/e^,t5} 


I ^ ”1” 1 11 II ^11 

r llL”([o,oo)xnx5i) — Ir llL“([o,oo)xnxSi) 


< 


Hence, we naturally have 

( 2 . 22 ) 


,k+l I 


\ _ g-(H-A) max{t/e^,tt} 

TTa 

1 n 


e"(i+^)(‘'>-'’)ds 


I ([0,oo) X 0x51 ) 


< 


I L°°([0,oo)xOx5i) — ^ II IIl°°{[0,oo)xOx51) 

Thus, this is a contraction sequence for A > 0. Considering = u\, we have 

fe-i 


(2.23) 


lL“([0.oo)xnx5i) 


< 


1 


1 + A 


A||Lcx=([0_oo)xnxSi) ’ 


for fc > 1. Therefore, converges strongly in L°° to a limit solution u\ satisfying 


(2.24) 


hAll 




< 


l + A 


L°°([0,oo)xnxSi) — / y 11^ llL°°([0,oo)xnx5i) — 
fe=l 


L“>([0,oo)xnx5i) ■ 


Since u\ can be rewritten along the characteristics as 
(2.25) u\(x,w) 

= - etbW,w)e~‘^^+^'>*'’ + J f(t-e'^(tb-s),x-e(tb-s)w,w)e~^^+^^‘^*‘’~‘''>ds^ 

+'i-{t<eH,} (h(x - (etw)/e^, u;)e-(i+^)‘/^' + ^ ' f(eh, x - e(t/e^ - s)w, zi;)e-(i+^)(‘/^'-")ds^ , 
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based on Lemma O we can directly estimate 

(2.26) lhA||Loo([0,oo)xnx5i) - llfflli“([0,oo)xr-) + ll^llL~(nx5i) + ll/llL°°([0,oo)xnx5i) • 

Combining (12.241) and (I2.26L we can easily deduce the lemma. □ 

2.2. Estimate. It is easy to see when A —>■ 0, the estimate in Lemma 12 .2 1 blows up. Hence, we need to 
show a uniform estimate of the solution to the penalized neutron transport equation (12.1611 . 

Lemma 2.3. (Green’s Identity) Assume f{t, x,w), g{t,x,w) £ L°°{[0,oo)xnxS^) and dtf+ 'w-y xf, dtg + 
w ■ 'S/xg G L^([0, oo) X H X 5^) with f,g€ L^([0, oo) x L). Then for almost all s,t £ [0, oo), 

(2.27) [ [[ ({dtf+ w ■Wxf)g + (dtg+ w ■Wxg)f)dxdwdr 

Js JJnxS^ V / 

= [ [ f9dldr+ ff f{t)g{t)dxdw - ff f{s)g{s)dMw, 

Js Jr JJqxS^ JJqxS^ 

where d'y = {w ■ n)ds on the boundary. 

Proof. See O Chapter 9] and [3]. □ 

Lemma 2.4. The solution u\ to the equation h2. 1 61) satisfies the uniform estimate in time interval [s,t], 

(2.28) 

II^A|lL2([s_j]xnx5i) - C'(^) ^ \\u\ - WA|lL2([s_t]xnxSi) + ll/llL2([s_t]xnx5i) + II^A|lL2([s_j]xr+) 

ll5llL2([s,t]xr-) ^ “ e^G'(s), 

where G(t) is a function satisfying 

(2.29) G{t) < C{n) ||1^A(^)||/^2(f2xSl) > 

for 0 < A << 1 and 0 < e << 1. 

Proof. We divide the proof into several steps: 


Step 1: 

Applying Lemma [27^ to the solution of the equation p.l6|l . Then for any (p £ L°“([0, oo) x H x 5^) satisfying 
edtp + w ■ yxfi £ -b^([0, oo) X H X 5^) and (p £ L^([0, oo) x L), we have 


(2.30) 


A 


nxsi 


u\(p - € 


nxsi 


dtpux - e 


0x51 


{w- S/x(p)u\ + 


0x51 


(UA - Ux)p 


f [ uxfpd'j - [ [ ux{t)p{t) + f f ux{s)p{s) + f ff fp. 

Js Jr ddOx5i JJQxS^ Js JJnxS^ 


Our goal is to choose a particular test function p. We first construct an auxiliary function ({t). Since 
ux{t) £ L°°{n X S^), it naturally implies ux{t) £ L°°(rt) which further leads to uxit) £ L^(H). We define 
C(t,x) on O satisfying 

AC(t) = uxft) in O, 

(ft) = 0 on dU. 

In the bounded domain O, based on the standard elliptic estimate, we have 
(2.32) ||C(t)|| C'(0)|fyA(t)|| L2(0) ■ 


(2.31) 


Step 2: 

Without loss of generality, we only prove the case with s = 0. We plug the test function 
(2.33) Pit) = -w-yxC{t) 

into the weak formulation (12.301) and estimate each term there. Naturally, we have 
(2-34) < C ||C(i)llpri(o) — C'(H) ||'UA(t)|li2(Q). 
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Easily we can decompose 


(2.35) - e 


0 JJQxS^ 


{w-S/x(t>)u\ = -e 


0 


{w ■ Vx4>)u\ - e 


0 JJQxS^ 


(w ■ Vx(l>)iu\ - Ux). 


We estimate the two term on the right-hand side of (12.351) separately. By ()2.31l) and (I2.33L we have 

(2.36) 


—e 


0 JJQxS^ 


{w-S/^(j))ux = e 


= e 


/ // Ux[wi{widiiC,+W2di2C,) +W2{widi2C,+W2d22(,) 

Jo JJnxS^ V 

UxiwldiiC + W2d22C 


0 JJQxS^ 


= eTT Ux{diiC + 0220 

Jo Jci 

= eTT l|wA|li2([o,t]xn) 

In the second equality, above cross terms vanish due to the symmetry of the integral over 5^. On the other 
hand, for the second term in (12.351) . Holder’s inequality and the elliptic estimate imply 

(2.37) -e ' '' ^ '• " - " ' ' 


Iff2(n) 


{w -W^Oiux-Ux) < C'(f^)e||wA-UA||i2(rotlxnx5i) ( 

Jo J JQxS^ \Jo 

< C{n)e ||ma - WA|lL2([o,t]xnx5i) II^A|li 2 ([o_t]xnxSi) ■ 

Based on ()2.32p . (12.341) . the boundary condition of the penalized neutron transport equation (j2.16|l . the trace 
theorem. Holder’s inequality and the elliptic estimate, we have 

(2.38) 


0 Jr 


ux^id'y = e / / UA^dy -I- e 


0 Jr+ 


UAd’dy 


0 ^r- 


< (7(0) l^e ||MA|li2([o_(]xr+) ll“^llL2([o,t]xnx5i) + ^ ll5llL2([o,t]xr-) ll^^llL2([o,t]xnx5i) 

(2.39) A r //' UA</> = xf [[ uxcl^ + xf [[ {ux - ux)^ = X f f f {ux-ux^ 

Jo JJqxs^ Jo J Jqxs^ Jo JJqxs^ Jo JJqxs^ 

< C{n)X ||uA|lL2([o,t]xnx5i) 11“^ “ '“A|lL2([o^t]xnxSi) I 
(2-40) f jj («-A - Ua)()> < < 7 ( 0 ) ||MA|li;,2([o_t]xnx5i) 11^'’' “ “'’'llL2([0,t]xnx5i) ’ 

Jo JJqxS^ 

(2-41) / // /<^ ^ <^(^) ll“A|lL2([o,t]xnxsi) ll/llL2([o,t]xnxsi) ■ 

Jo JJClxS^ 

Note that we will take 

(2.42) -e^ [[ uximt) +0 [[ uxiom =o(G{t) - G(0) 

JJqxs^ JJqxs^ \ 

where G{t) = — u\{t)4>i't)- Then the only remaining term is 

J JqxS^ 


(2.43) 


dt(l)ux = -e 


f f f dt(j){ux - Ux) 

Js JJqxS^ 


0 J J JQ,xS^ 

< ll^tVC||_f^2([o_t]xnxsi) “ '“A|lL2([o_t]xnx5i) • 


Now we have to tackle ||9tVC||i2([o,t]xnxSi 
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Step 3: 

For test function 4’{x, w) which is independent of time t, in time interval \t — S, t\ the weak formulation in 
(I2.30|) can be simplified as 

(2.44) \[ ff ux(j)-e[ [[ {w-Vx4>)u\+ [ [[ {u\-ux)4> 

Jt-sJJnxS^ Jt-sJjQxS^ Jt-sJJnxs^ 

ux4>d'y - [[ ux{t)4> + [[ uxit - 6)4> + [ [[ f4>. 

jJnxS^ jJnxS^ Jt-sJJnxS^ 


= -^LsJr 

Taking difference quotient as <5 —>■ 0, we know 


(2.45) 

Then (12.441) can be simplified into 

(2.46) [[ dtux{t)(j) 

J JQxS^ 


[[ ux{t)(j) - [[ ux{t-5)4> 

JJnxs^ JJqxs^ 


[[ dtux{t)(j). 

J Jnxs^ 


-A /'/' ux{t)(j) + e[[ (w ■ Vx4>)ux{t) - (( {uxit) - ux{t))(j) 

JJnxs^ jJnxs^ jJnxs^ 

-e [ ux{t)(l)d'y+ ff f{t)4>. 

Jr JJnxS^ 


For fixed t, taking (p = $(a;) which satisfies 

r A$ = dtux{t) in n, 

^ ’ \ $ = 0 on 

which further implies $ = dtC- Then the left-hand side of (I2.46|) is actually 


(2.48) 


LHS = 


ff ^dtux{t) = ff ^dtux 

JJqxs^ JJqxs^ 

ff |V4>|' 

J Jnxs^ 




nx5i 


— ll^tVC(t)||i2(f2xSi) • 

By a similar argument as in Step 2 and the Poincare inequality, the right-hand side of (j2.46l) can be bounded 


as 


(2.49) 

RHS < ||5tVC(t)||i2(f2x5i) ^ l|■^A(0 ~ '^A(t)|li2(Qx5i) + ll'^'>'(0llL2(nxSi) + ll/(0llL2(nxSi) ^ ■ 
Therefore, we have 

(2.50) ||9tVC(t)|| i2(f2x5i) ^ ll'^>(0 ~ '^A(t)||j;,2(f2x5i) + ll'^'*'(^)llL2(nx5i) + II/(^) 1^2(0x51) ' 

For all t, we can further integrate over [0, t] to obtain 
(2-51) ll^t^C(^)llL2([o,t]xnx5i) 

< \\ux{t) - 'WA(t)|lL2([o_t]xnx5i) + ^ II^A(t)|lL2([o,t]xnx5i) + 11/(0 llL2([o,t] xHxSi) • 


Step 4: 

Collecting terms in (12.361) . (I2.37p . (12.381) . (I2.39L (12.401) . (I2.41L (12.421) . (12.431) . and (I2.51L we obtain 
(2.52) e ||MAl02([o_i]xnx5i) 

< C{il.) ^(1 -I- e -b A) ||ma - '«A|lL2([o_t]xnx5i) + ^ l|■“A|lL2([o_t]x^+) + ll/llL2([o,t]xnx5i) + ^ ll5llL2([o_t]xr-) 
+e‘^Git) - e‘^G{0). 
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When 0 < A < 1 and 0 < e < 1, we get the desired uniform estimate with respect to A. 


□ 


Theorem 2.5. Assume f{t,x,w) S L^([0,oo) x x 5^), h{x,w) S x S^) and e^'^*g{t,xo,w) S 

L^([0,oo) X F”) for some Aq > 0. Then for the unsteady neutron transport equation there exists Ag 

satisfying 0 < Ag < Ag and a unique solution u{t,x,w) € L^([0, oo) x fix 5^) satisfying 


(2.53) 


-.At ^ 


-.At. 


pl/2 


'“'lL 2 ([o_oo)xr+) ■*" Ir Ir 


L2([0,oo)xnxSi) 


< (7(12)^^2 116 /|L2([g_oo)xnx5i) Il^lli=(nx5i) + gi/2 11® 5||^2([o_3c)xr-) ^ ’ 


for any 0 < A < Ag. When Ao = 0, we have Ag = 0. 

Proof. We divide the proof into several steps: 

Step 1: Weak formulation. 

In the weak formulation (12.301) . we may take the test function (jj = u\ to get the energy estimate 


(2.54) 




Je^l|wA(t)|li2(ox5i) - 


- i|2 


lL2(nx5i) + ll“A - '«A|li2([g_t]xnxSl) 


III 

Jo JJQ: 


fu\. 


10 J JQxS^ 
Hence, this naturally implies 

< 


(2.55) 


6 ll“A|ll2([o_t]xr+) + “6” II^aI 


1 


lL2(nxSi) + ll“A - ''^A|li2([o t]xnx5i) 


/wA + ^e ll^llL=(nxsi) 


ri 


/o JJQxS^ 2 

On the other hand, we can square on both sides of (12.281) to obtain 
(2.56) 6^ II^A|li^2([o tjxOxSi) 

|2 , II J-||2 , 2 II. I|2 


L2([o,t]xr-) ■ 


< (7(0)( ||ua ^A|li2QQ (jxnxsi) + ll/llL2([o,t]xnx5i) "*■ 6 l|wA|li2([Q_(]xr+) + ^ ll5llL2([o,t]xr-) 


+e^ ||uA(t)||2,2(QxSi) 


Il2(OxS1) 


Multiplying a sufficiently small constant on both sides of (I2.56P and adding it to (12.551) to absorb ||ua|1^2(p+), 

l|wA(t)|li 2 (oxsi) and ||uA - WA|li 2 (nxsi): ^e deduce 

2 2^22 2 
(2.57) e ||MA|lL2([g t]xr+) + 6 I^a( 011^2(0x51) + ^ l|wA|lL2([o_t]xnx5i) + II^a - ■aA|lL2([o tjxnxsi) 

t 


< cm ii/iii2([g,]xnx5i) +IJ /«A 


lL2(nx5i) 


lL2([o,t]xr-) ■ 


Hence, we have 


(2.58) 


II^A|li2([o^t]xr+) + II^a(0IIl2(Qx 51) + ll“A||p2([g_t]xnx51) 


< C'(H)( ||/||p2([o^t]xnx5i) + 


10 JJqxS^ 

A simple application of Cauchy’s inequality leads to 

1 


fu\ + e ||/i|lz,2(nx5i) + 6 ll5llL2([g_t]xr-) )• 


(2.59) 


fux < 

to JJQxS^ 4Ce 


L2([0,t]xnx5i) + Ce l|■aA|lL2([o_(]xnx5l) ■ 
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Taking C sufficiently small, we can divide (I2.58|) by to obtain 


(2.60) 


- II^A|lL2([o_t]xr+) + II^A(i)|li2(f2x5i) + l|■“>ll 


< 


lL2([o,t]xnx5i) 


lL2(nx5i) 


L2([o,t]xnx5i) 

1 ,, 1,2 

e 


lL2([o,t]xr“) 


Step 2: Convergence. 

Since above estimate does not depend on A, it gives a uniform estimate for the penalized neutron transport 
equation (12.161) . Thus, we can extract a weakly convergent subsequence ua —t u as A —?> 0. The weak lower 
semi-continuity of norms IMIl 2 ([o t]xnx 5 i) l|■|lL 2 ([o t]xr+) implies u also satisfies the estimate (12.601) . 
Hence, in the weak formulation (I2)30p . we can take A —>■ 0 to deduce that u satishes equation (12.ip . Also 
u\ — u satisfies the equation 

(2.61) 


^dtiux -u) + ew ■ Vx{u\ -u) + {u\ - u) - {u\ - u) = -Xu\ 




(ux — u)(0,x,w) = 0 in H, 

{ux — u){xo,w) = 0 for xo € dfl a,ndw-n<0. 

By a similar argument as above, we can achieve 


(2.62) 


\\ux - u||i 2 ([o,t]xnxsi) ^ C{fl) ( l|MA||i2([o_t]xnx5i) ] • 


A 


When A —>■ 0, the right-hand side approaches zero, which implies the convergence is actually in the strong 
sense. The uniqueness easily follows from the energy estimates. 


Step 3: Decay. 

Let V = e^*u. Then v satisfies the equation 

{ e^dtv -I- ew ■ VxV + v — v = f + Ae^u in D, 
v{0,x,w) = h{x,w) in Q, 

v{xo,w) = e^*g{t,Xo,w) for Xq & dV, a.ndw-n<0. 
Similar to the argument in Step 1, we can obtain 

(2.64) 1 Hull 

L2([o,t]xr+) + Mm 

L2(nx5i) + ll^llL2([o,t]xnx5i) 


< C(D)(l||e"V"' 


1 


_l_ \z 4 

.4 ir J llL2(ro #1x0x51) ^ ^1 


n - ,lL2([o_t]xOx5i) 
Then when A is sufficiently small, we have 


L2([o,t]xOx5i) 


llL2(nx5i) + g II® f^lL2([o_t]xr-) 


(2.65) 


1 2 

~ ll^llL2([o,t]xr+) ^ 


II'*^(0I1l2(Qx51 

L2([o,i]xOx5i) 


Ikll 


L2([o,t]xOx5i) 


1 


Xt, 


+ ll“llL2(nx5i) + II® 5|lL2([o,t]xr-) )’ 


which implies exponential decay of u. 


2.3. L°° Estimate. 


□ 


Theorem 2.6. Assume f(t,x,w) G ^°“([0:Oo) x H x 5^), h{x,w) S x 5^) and e^°*g{t,xo,w) S 

L°“([0, oo) X T”) for some Aq > 0. Then for the unsteady neutron transport equation there exists Ag 

satisfying 0 < Ag < Aq and a unique solution u{t,x,w) € L°°([0,oo) x fl x 5^) satisfying 

( 2 . 66 ) 


® ^IL“([o,oo)xr+) + II® ^(^)IL°°(nx5i) II® '“IL“([o,oo)xnx5i 




At 


< c(fi)(4ll®"7 


I L^([0,oo)xr2XtS 




2 ll^llL2(nx5i) + 5/2 ll®”'’5|lL2(io_^)xr-) 


At ^ 


+ ||®^7 


I L°° {[0,oo)xQxS^) 


lL«^(nx5i) 


+ l|e^7| 


L“([0,oo)xr-) /> 
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for any 0 < A < Ag. When Aq = 0, we have Ag = 0. 

Proof. We divide the proof into several steps to bootstrap an solution to an L°° solution: 

Step 1: Double Duhamel iterations. 

The characteristics of the equation (12. ip is given by (12. lip . Hence, we can rewrite the equation (12.11) along 
the characteristics as 

(2.67) u(t,x,w) 

= - etbW,w)e~^’’+ J (u + f){t - e^(tb - s),x - e(tb - 


+ l{t<e 


^h{x — (etw)/, w)e + J {u + f){e‘^s,x — e{t/e'^ — s)w,w)e 


where the backward exit time tb is defined as (12.141) . For the convenience of analysis, we transform it into a 
simpler form 


( 2 . 68 ) 


u{t, X, w) 


= ^{t>eHb}9{t-e‘^tb,x-etbW,w)e *” + - {<^tw)/e'^,w) 


2 


rt/P 


+ 


1 r*/ 


f{e^s,x — e{t/e^ — s)w,w)e '’Ids 


f ( f u{e^s,x — e{t/e'^ — s)w,wt)dwi\e d/” '’Ids. 

Jt/e^ — th^it/e^) \J — -K / 


Here a A 6 denotes min{a, &}. Note we have replaced u by the integral of u over the dummy velocity variable 
Wt- For the last term in this formulation, we apply the Duhamel’s principle again to u{e'^s, x—e{t/e'^ — s)w, Wt) 
and obtain 

(2.69) u{t,x,w) 

= - e‘^tb,x - etbW,w)e~*>’ + - {etw)/ 

f{e^s,x — e{t/e^ — s)w,w)e~^*^'^ 


ft/e 


h/e^-ti,A{t/e'^) 

PIP 


-- [ 
2711 


t/P 


/ 1! — tbAit j 

1 ^ 


27r 


H j — tbAit! 


J '^{s>sb}g{e^{s - Sb),x- e{t/e^ - s)w, wt)e~'''’dwt'^ “'*^ds 

J '^{s>sb}Hx- - s)w- eswt,u;t)e“Mwt^e“d/'’''-®)ds 

( f [ f{e'^r,x — e{t/e^ — s)w — e{s — r)wt,wt)e~^^~'^^drdwtje~^*^'^~^^ds 

\J—7TJs — SbAS / 


L /•"“ 

271 Jt/e^-ttA(t/e^) 


( [ f u{e^r,x — e{t/e^ — s)w — e{s — r)wt)e ’’^drdict'je ®^ds. 

\ J —7T J S — SfjAS J 

where we introduce another dummy velocity variable Ws and 

(2.70) Sb{x,'w,s,wt) = inf{r > 0 : (a? — e{t/e^ — s)'w — erwtTWt) G F”}. 


Step 2: Estimates of all but the last term in (I2.69F 
We can directly estimate as follows: 

(2.71) ehb,x-etbw,w)e~^'’\ < ||g|lLoo([o^t]>,r-) > 
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(2.72) 

(2.73) 

(2.74) 

(2.75) 

(2.76) 

< 


11 — j €■ 
tie 


- {etw)/e^,w)e 
f(e‘^s,x — — s)w,w)e~^*^'^ 


^ ll^llL~(nx5i) > 


^ ll/llL“([0,t]xnxSi) ’ 


^ [ ( [ l{s>Sb}9i<^'^is - Sb),x - eit/e"^ - s)w,Wt)e ®Ms 

Jt/e^-ttA(t/e2) \j--rr J 


< 


L°°([o.t]xr-) ’ 
1 7*/*^ 


^ ( ( f l{s>Sb}Hx - - s)w - eswt,Wt)e 

Jt/e^-tbA(t/e^) \J-n ) 


< 


tf e'^ — tbA{tje 

lL“(nx5i) ’ 


^ [ ( f f f{e^r,x — e{t/e^ — s)w — e{s — r)wt,wt)e ’'MrdiPt^e 

J tje^ — tbAitje^) \J—7rJs—SbAs J 


lL°®([o,t]xnx<si) • 


Step 3: Estimates of the last term in (12.691) . 

We can first transform the last term I in (12.691) into 


(2.77) 

|/| < 


f iff \u{e^r,x — e{t/e^ — s)w — e{s — r)wt)\e '’MrdiUt^e '’Ms 

^lle^ — thAitje^') \J—7rJs—SbAs / 


< 


1 Z*^/ 

27r 


t/e^ —tbA(t/e^) 
tb / /*7r r^b 


0 \ J —t: ^0 


\u{e^{r* + s* + t/e'^ — tb — Sb),x — e{tb — s*)w — e{sb — r*)wt)\e ’’ Mr*dw;t^e '’Ms*, 

by substitution s —>■ s* = (s — t/e^ + tb) and r —>• r* = (r — s + Sb)- Now we decompose the right-hand side 
in (12.771) as 


f'tb I* P^b P^h 


rtb 


(2.78) ///=/// +/ / / =h+h, 

Jo Js^Jo Jo Js^Jsi,-r*<S Jo Js^Jsi,-r*>S 

for some ^ > 0. We can estimate Ii directly as 

f'tb / P^b 

\L°°{[0,t]xflxS 


(2.79) h < f e"(‘'>"’'*)Z f 
Jo \ Jn 


Ullroo^rn twos. ci', dr* ds* < ^ llul 


max{0,Sb —(5} 


L°°{[0,t]xQxS^) • 


Then we can bound I 2 as 
(2.80) 

h < C 


tb p ^max{0,S5 —(5} 


^0 JS^Jo 

|'«(e^(r* -I- s* +tfe^ — h — Sb),x — e{tb — s*)w — e{sb — r*)wt)\ e“^'’’’“’' )(;Jr*d'u;jds*. 

By the definition of tf, and Sf,, we always have e^(r*-|-s*-|-t/e^—tt,—Sb) £ [0, t] and a:—e(tf,—s*)ii;—e(s&—r*)ii;t £ 
n. Hence, we may interchange the order of integration and apply Holder’s inequality to obtain 

(2.81) 
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rtb 


max{0,S5 —(5} 


h < C I I I + s*+t/e'^-tb-Sb),x-e{tb-s*)w - e{sb-r*)wt)e ’’*)■ 

Jo Js^ Jo 

|w(e^(r* + s* + t/— tb — Sb), X — e{tb — s*)w — e{sb — r*)'Wt) \ dr*dwtds* 

tb 


< C 


IQ JS^ JQ 
(•tb 


max{0,Sb —5} \ 1/2 

e-(^'>-^'‘)-d»-^’')dr*dwtds* 

max{0,S5 —(5} 


/O JS^ Jo 
^2/_* 


+ s*+t/e'^ - tb-Sb),x - e{tb-s*)w - €{sb-r*)wt)e 

\ 1/2 

|M^(e^(r* + s* + //e^ — tb — Sb),x — e{tb — s*)w — €{sb — r*)wt) \ dr*dwtds* j 

Note Wt G 5^, which is essentially a one-dimensional variable. Thus, we may write it in a new variable ip as 
Wt = {cosil),smip). Then we define the change of variable [0,t/e^] x [—7r,7r) x R —>• [0, t] x : {s*,il),r*) —>■ 
{t',yi,y2) = (£^(r* -b s* -I- t/e^ - tb - Sb),x- e(tb - s*)w - e(sb - r*)wt), i.e. 

( t' = e'^{r* + s* +t/e'^ - tb - Sb), 

(2.82) } = xi - e{tb - s*)wi - e{sb - r*) cosip, 

{ y 2 = X 2 - e{tb - s*)w 2 - e(sb - r*) s'mip. 

Therefore, for Sb — r* > 6, we can directly compute the Jacobian 

(2.83) 

^ 0 e2 

ewi e{sb — r*) sin Ip e cos ip 
€W 2 —e{sb — r*) cosip esinip 

Thus, in order to guarantee the Jacobian is strictly positive, we may further decompose I 2 into / 2 ,i + 12,2 
where in / 2 ,i, we have wi cosip + 1 x 2 sin'i/i > 1 — J and in J 2 , 2 , we have wi cosip + 1 x 2 sin'i/i < 1 — J. Since 
w = {wi,W 2 ) G 5^, based on trigonometric identity, we obtain 

(2.84) 


d{t',yuy2) 


d{s*,ip,r*) 



= e'^(sb — r*) 1— (wi cos ip + 1 x 2 sin ip) j. 


I2,i < C 


^max{0,Sb —(5} 


'0 Jw-wt>i^—S J0 

|M(e^(r* -b s* +tle'^ — h — Sb),x — e{tb — s*)w — e(sb — r*)r(;t)| )-(tb-s )(^r*dixtds* 

< <5 ||M|lLoo([o_t]xnx5i) ■ 

Hence, we may simplify (12.811) as 

l-t I- I- ^ \l/2 

(2.85) /2,2 < ' 




Then we may further utilize estimate of u in Theorem 12.51 to obtain 

C 

^^2,2 < ^ ll'“llL2([o,(]xnx5i) 


( 2 . 86 ) 


< 


g ll•/'llL 2 ([o,t]x^xSl) + g 2 ll^llL 2 (nxSi) + g 5/2 ll5llL2([0,t] xT") ^ ' 


Step 4: L°° estimate. 

In summary, collecting (I2.71L (j2.72l) . (12.731) . (j2.74l) . (12.751) . (I2.76L (I2.79P and (12.861) . for fixed 0 < J < 1, we 
have 

(2.87) \u{t,x,w)\ 

< S ||u|lLoo([o_t]xnxSi) ”1 ^ f ^ ll/llz, 2 ([o,t]xaxSi) + ^ II^IIl2(Ox51) + ^572 ll 5 'llL 2 ([ 0 ,t]xr-) 


+ ( ll/llL°°([0,t]xax5i) + ll^llL°°(nx5i) + ll 5 llL“([ 0 ,t]xr- 


(tb-s*) 
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Then we may take 0<(5<l/2to obtain 

(2.88) \u(t,x,w)\ 

1 „ „ _/ 1 


— 2 ll^lli“([o,t]xnxsi) + C(^^)( 4 ll/llL2([o,t]xnxsi) + 2 ll^llT2(nxsi) + 5/2 ll 3 llL 2 ([o,t]xr-) 


+ ^ ll/llL~([0,t]xnx5i) + ll^llL“(nx5i) + llffllL“([0.t]xr-) 
Taking supremum of u over all (t, x, w), we have 
(2.89) IKIlL“([o,t]xax5Q 


— 2 ll^lli“([o,t]xnxsi) + C'(r2)l ll/llL2([o,t]xnxsi) + 2 ll^llL2(nxsi) + 5/2 ll 5 llL 2 ([o,t]xr-) 


+ ll/llL=°([0,t]xnx5i) + ll^llL°°(nx5i) + llffllL“([ 0 .t]xr-) 

Finally, absorbing ||w|lLoo([o,t]xnxSi)^ get 

(2.90) ll^llL~([o,t]xnxSi) — ll/llL2([o,t]xax5Q + ^ 11^11^2(0x51) + llffllL2([o,t]xr-) 

+ ll/llL“([0,t]xnx5i) + II^IIl“(Ox 51) + ll5llLoo([o_t]xr- 


) ■ 


Step 5: L°° Decay. 

Let V = e^‘u. Then v satisfies the equation 

{ e^dtv + ew ■ VxV + (1 — Ae^)n — v = / in D, 

v(0,x,w) = h(x,w) in fl, 

v{xo,w) = e^*g{t,xo,w) for xq G dil andw-n<0. 

By a similar argument as in Step 3 and Step 4, combined with the decay, we can finally show the desired 
estimate. □ 


2.4. Maximum Principle. 

Theorem 2.7. When / = 0, the solution u(t,x,w) to the unsteady neutron transport equation satisfies 
the maximum principle, i.e. 

(2.92) mi{g{t,xo,w),h{x,w)} < u{t,x,w) < sup{g(t,xo,w),h{x,w)}. 

Proof. We claim that it suffices to show u{t,x,w) < 0 whenever g(t,xo,w) < 0 and h{x,w) < 0. Suppose 
the claim is justified. Then define 

(2.93) m = mi{g{t,xo,w),h{x,w)}, 

(2.94) M = sup{g{t,xo,w),h{x,'w)}. 

We have ui = u — M satisfies the equation 

{ e^dtUi + ew ■VxUi + ui — ui = 0 in [0,oo)xD, 

ui{0,x,w) = h{x,w) — M in D 

ui{t,xo,w) = g{t,xo,w) — M for w-n<0 and Xg G dfl, 

Hence, h — M < 0 and g — M < 0 implies ui < 0, which is actually u < M. Similarly, we have U 2 = m — u 
satisfies the equation 

{ e^dtU 2 + ew ■ VxU 2 + U 2 — U 2 = 0 in [0, 00 ) x D, 

U 2 { 0 ,x,w) = m — h{x,w) in D 

U 2 {t,xo,w) = m — g{t,xo,w) for w • n < 0 and xq G dQ, 

Hence, m—h<0 and m — g < 0 implies U 2 < 0, which is actually u> m. Therefore, the maximum principle 
is established. 
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We now prove the claim that if g{t,xo,w) < 0 and h{x,w) < 0, we have u{t,x,w) < 0. We first consider 
the penalized neutron transport equation 

{ Xu\ + e^dtUx + ew ■ VxUx + u\ — u\ = 0 in [0,oo)xn, 

ux{0,x,w) = h{x,w) in 

ux{t,XQ^w) = g(t,xo,w) for Xq G dfl a.ndw-n<0. 

with A > 0. Based on Lemma [221 there exists a solution ux{t,x,w) G L°°([0, oo) x n x 5^). We use the 
notation in the proof of Lemma |2.21 Define an approximating sequence where = 0 and 

( Xu’l + e'^dtu'l + ew ■'Vxu'l + u'l — = 0 in [0, oo)xD, 

(2.98) < u*(0,af, iZl) = h{x,w) in D 

t u^{t,xo,w) = g{t,xo,w) for xq G dV, andw-n<0. 

By Lemma [^m this sequence is well-defined and ||^^a||loo([o oo)xnx 5 i) ^ rewrite equation (j2.18l) 

along the characteristics as 

(2.99) 


u\{t,x,w) 

= l{t>e2t,} - eHb, X - etbW, u;)e"(^+^^*'> 4 

+I{t<e2t,} ^/i(x - (etw)/e^iZ;)e"(^+^)‘/'^' -h 






-M- 


{t - e^{tb - s),x- 
^(e^s, X — e(t/e^ — 


0 


where 


t{tb-s)w, iZ;)e-(i+^)(*'>-")ds 

s)ui,u;)e-(i+^)(‘/^'-")dsY 


(2.100) tb{x,w) = inf{s > 0 : {x — €sw,w) S L }. 

Since u\{t,x,w) < 0 naturally implies u\{t,x) < 0, we naturally have u'l^tjXjW) < 0 when g{t,xo,w) < 0 
and h{x, w) < 0. In the proof of Lemma [2.21 we have shown u\ —>■ ux in L°° as /c —>■ oo. Therefore, we have 
uxit, X, w) < 0. Based on the proof of Lemma [2.51 we know ua —>■ w in as A —>■ 0, where u is the solution 
of the equation (EH). Then we naturally obtain u < 0. Also, this is the unique solution to the equation 
(EH). This justifies the claim and completes the proof. □ 


Theorem 12.71 naturally leads to the L°° estimate of the equation (12.11) . 

Corollary 2.8. Assume h{x,w) G L°°{V, x 5^) and g{t,xo,w) G L°°{[0,oo) x T”). Then for the unsteady 
neutron transport equation \2.1]) with / = 0, there exists a unique solution u{t,x,w) G L°°([0, oo) x D x 5^) 
satisfying 

(2.101) ||M|lL=o([ 0 _oo)xr+) + lk(0llL~(nxSi) + ll^llL“([0.oo)xnx5i) - ll^llL“(nx5i) + ll5llL=o([0,oo)xr-) • 

2.5. Well-posedness of Transport Equation. Combining the results in Theorem 12.61 and Corollary EH 
we can show an improved L°° estimate of the equation EH- 

Theorem 2.9. Assume f{t, x, w) G L°°([0, oo)xDx5^), h{x, w) G L°°(Dx5^) and g{t, xo,w) G L°°([0, oo)x 
r“). Then for the unsteady neutron transport equation there exists a unique solution u{t, x, w) G 

L°“([0, oo) X fl X 5^) satisfying 

(2.102) ll^llL~([o,oo)xr+) + ll^(0llL~(nxSi) + ll^llL“([o,oo)xnx5i) 

< (7(11)^^ ll/llL2([0_oo)xnx5i) + ll/llL“([0,oo)xnx5i) ^ ll^llL”(nx5i) + llffllL“([ 0 ,oo)xr-) ■ 

Proof. Since the equation EH is a linear equation, then we can utilize the superposition property, i.e. we 
can separate the solution u = ui -I- M 2 where ui satisfies the equation 

f e'^dtUi + ew ■ X/xUi + ui — ui = 0 in [0, oo)xf2, 

ui{0,x,w) = h{x,w) in D 

ui(t,XQ,w) = g{t,xo,w) for w-n<0 and xq G dfl, 


(2.103) 
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and U 2 satisfies the equation 

{ e'^dtU 2 + ew ■'VxU 2 + U 2 — U 2 = f(t,x,w) in [0,oo)xr2, 

U 2 { 0 ,x,w) = 0 in 

U 2 {t,Xo,w) = 0 toi w ■ n < 0 and Xq & dV,, 

Note that the data in (I2.103|) and (I2.104|) satisfy the compatibility condition (11.31) . Therefore, we can apply 
the previous results in this section. Corollary 12.81 yields 

(2.105) hill^ 

“([ 0 ,oo)xr+) + ll'“l(0llL“(nx5i) + ll^lllL“([0,oo)xnx5i) — ll^llL“(nx5i) ll5llL°°([0,oo)xr-) • 

Also, Theorem 12.61 leads to 

ll^2|lioo([Q_3c)xr+) + ll''^2(0llL°°(nx5T ll^2|lioo([Q_oo)xnx5Q 

< C{^) ll/llL2([0,oo)xnx5i) + ll/llL“([0,oo)xnx5i) ^ • 

Combining (I2.105P and (12.1061) . we have the desired result. □ 


Finally, we can apply Theorem 12.91 to the equation and obtain Theorem o 

Theorem 2.10. Assume g{t,XQ,w) € L°°([0,oo) x T”) and h{x,w) € L°°(fl x S^). Then for the unsteady 
neutron transport equation urn, there exists a unique solution u^{t,x,w) € L°“([0,oo) x fl x 5^) satisfying 

(2-107) l|■“'^llL“([o.oo)xnx5l) — ll^llL“(nxSi) + ll5llL~([o.oo)xr-) • 


3. Asymptotic Analysis 

In this section, we construct the asymptotic expansion of the equation dm. 

3.1. Discussion of Compatibility Condition. The initial and boundary data satisfy the compatibility 
condition 

(3.1) h(xo,w) = g{0,xo,w) for w ■ n < 0. 

Then in the half-space w • n < 0 at {0,Xq,w), the equation 

(3.2) e^dtu^ + ew ■ SIxu'^ -\-u'^ — = 0, 

is valid, which implies 

(3.3) e^dtg{Q,xo,w) + ew-y xh{xo,'w)-\-h{xQ,w)— h{xo) = 0. 

In order to show the diffusive limit, the condition (13.31) holds for arbitrary e. Since g and h are all independent 
of e, we must have for ifi • fi < 0, 

(3.4) dtgi0,XQ,w) = 0, 

(3.5) w ■Vxh{xo,w) = 0, 

(3.6) h{xo,w) — h{xo) = 0. 

The relation (13.61) implies the improved compatibility condition 

(3.7) h{xQ,w) = g{0,XQ,w) = Cq for w ■ n < 0, 

for some constant Cq. This fact is of great importance in the following analysis. 
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3.2. Interior Expansion. We define the interior expansion as follows: 

oo 

(3.8) W{t,x,w) 

k=0 

where can be defined by comparing the order of e via plugging (13.81) into the equation (11.11) . Thus, we 
have 


(3.9) 

(3.10) 

(3.11) 


U^-Uo = 0, 

ut-ut = -^z;•v,c/o^ 

U^-UI = -w-V,Ul-dtU^, 


(3.12) 


U^k-UI = -w-yM_,-dtUU. 


The following analysis reveals the equation satisfied by {7|: 

Plugging (13.91) into (13.101) . we obtain 

(3.13) C/i" = C7^-'uJ-V,,C7o". 

Plugging (13.131) into (13.lip , we get 

(3.14) Ui - C7| + dtU^ = -w ■ V,(C7^ - w ■ = -w ■ + \w\'^ + 2wiW2d,,,,U^. 

Integrating (13.141) over w G S^, we achieve the final form 

(3.15) dtU^ - = 0, 

which further implies UQ{t,x,w) satisfies the equation 


(3.16) 


US 

dtUs - A,us 


US, 

0 . 


Similarly, we can derive Ul{t,x,w) satisfies 


(3.17) 


ui 

dtUl - A,Ul 


and U^{t,x,w) for fc > 2 satishes 


Ul-w-V^US, 

0 , 


(3.18) 


= U^,-w-V,UU-dtU^k-2, 


Uk 

dtUl-A,U^^ = 0. 

Note that in order to determine [/^, we need to define the initial data and boundary data. 


3.3. Initial Layer Expansion. In order to determine the initial condition for U^, we need to define the 
initial layer expansion. Hence, we need a substitution: 


Temporal Substitution: 

We define the stretched variable r by making the scaling transform for u'^(t) —>■ w'^(t) with r G [0,oo) as 


(3.19) 

which implies 


(3.20) 


du^ 


1 du^ 
dr 


In this new variable, equation (ini) can be rewritten as 


dru'^ + ew ■ VxU^ + = 0, 

u'^(0,x,w) = h(x,w), 
u'^(T,xo,w)=g(T,xo,w) for w ■ n < 0. 


(3.21) 
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We define the initial layer expansion as follows: 

OO 

(3.22) %%T,x,w) - '^e^%^k{T,x,w), 

k=0 

where can be determined by comparing the order of e via plugging (13.221) into the equation (13.211) . 
Thus, we have 


(3.23) 

^t'^1,0 + ^/,0 

-^to = 

0, 

(3.24) 

dr%\^ + 

-Kl = 

-W • Va;'^/0, 

(3.25) 

dr^j2 + ‘^1,2 

-%\2 = 

-w ■ 

(3.26) 


~ ‘^I,k = 



The following analysis reveals the equation satisfied by 
Integrate (|3.23|) over w £ we have 

(3.27) dr%^,o = 0. 


which further implies 


(3.28) 

Therefore, from (13.231) . we can deduce 

(3.29) %^o{t,x,w) 


%\e(T,x) = ‘^lo{0,x). 

= e“'^'^/o(0,af, iZ;) + f '^/q(s, a:)e®“'^ds 

Jo 

= e“^'^j(o(0,a?,u;) + (1 - e“'^)‘^/o(0, f). 


This means we have 
(3.30) 


dr^to 

'^/o(t, X, w) 


0 , 

e“'^'^j‘;Q(0,f, w) + (1 - e“'^)‘^j*;Q(0,x). 


Similarly, we can derive ‘Wff.{T,x,w) for fc > 1 satisfies 


(3.31) 




w ■ Vx^i^k-i ]dw 


‘^f ki'^,x,w) = e '^%'^^.{0,x,w) + J )(s,f,u;)e® ■"ds. 


3.4. Boundary Layer Expansion with Geometric Correction. In order to determine the boundary 
condition for C/^, we need to define the boundary layer expansion. Hence, we need several substitutions: 


Spacial Substitution 1: 

We consider the substitution into quasi-polar coordinates u’^{xi,X 2 ) —>■ with {p,9) £ [0,1) x [—7r,7r) 

defined as 


(3.32) 


f X 1 = (l — fj,)cos0, 
[ X 2 = (1 — /i) sin0. 


Here /r denotes the distance to the boundary dfl and 0 is the space angular variable. In these new variables, 
equation (HU can be rewritten as 


(3.33) 

f 29 u‘' ( . „\du‘' e f e 1 /" e 1 ^ n 

e —-e wi cos tl + ?ii 2 sm cZ —- wi sm 0 — W 2 cos 0 —— +u -/ udw = 0, 

\ at \ J I - ^Ji\ ) 00 2Tr Jsi 

I u''{0,fi,9,Wi,W2) = h{fl,0,Wi,W2), 

u'^{t,0,9,Wi,W2) = g{t, 0 ,Wi,W 2 ) for ici cos0 + u >2 sin0 < 0. 


Spacial Substitution 2: 
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We further define the stretched variable rj by making the scaling transform for —>■ u^(r],6) with 

( 77 , 0 ) S [0,1/e) X [—TT, tt) as 


(3.34) 

which implies 

(3.35) 

Then equation (HH) is transformed into 


[0 = 0 , 

du'^ 1 du'^ 
e dr] 


(3.36) 

{ 2 f n • ^ ( • /I 

e - wi cos 0 + W 2 sin 0 --- wi sin 0 — W 2 cos 

ot \ y ary 1 — 677 \ 

m''( 0, 77, 0, iZ;) = /i(77,0, wi,W2), 

u'^{t, 0, 9, wi, W 2 ) = g{t, 9,11)1, W 2 ) for wi COS0 + W 2 sin0 < 0. 


du'^ 


u — 


2ti 


u'^dw = 0 , 


Spacial Substitution 3: 

Define the velocity substitution for u'^ { 101 , 102 ) — t with ^ G [—tt, tt) as 

(3.37) I 

' ' ( 7C2 = —cos 4 . 

Here 4 denotes the velocity angular variable. We have the succinct form for (11.11) as 


, dii 


du' 


« ^+sin(^ + C)^- 


du’^ 


1 


(3.38) 


dt 




u''(O,77,0,4) = h{r],9,C), 

7i*^(t, 0,0,4) = g(t, 0,4) for sin(0 + 4)>O- 


Spacial Substitution 4: 

We make the rotation substitution for m'^( 4) —t u^{(j)) with (j) G [— 7 r, 7 r) as 
(3.39) (j) = 0 + 4, 

and transform the equation cu into 


, du’^ 


e~ + sin (f) 


du^ 


■ cost/ 


(3.40) 


du^ du'^ 
d4> d9 


u — 


2tt 


dt drj 1 — er] 

m"(O,? 7 , 0 ,(/) = h{r],9,(j)), 
u'^{t,0,9,4>) = g{t,9,4>) for sin^i > 0. 

We define the boundary layer expansion with geometric correction as follows: 


u'^dcj) = 0 , 


(3.41) 


‘^^{t, 77 , 0 , (/) - ^ e’^‘^^j,{t, 77 , 0 , (/), 


k=0 


where can be determined by comparing the order of e via plugging (|3.41|) into the equation (13.401) . 

Following the idea in m, in a neighborhood of the boundary, we require 


smt 


smt 


(3.42) 

(3.43) 

(3.44) sint/ 

(3.45) sinyi 


d% 


B,0 


d% 


dg 

1 

- e77 

dn,i 


e 

dg 

1 

- £0 

2 


e 

dg 

1 

- £0 

^^B,k 


e 


cost 


cos I 


B,0 


d(j) 


+ n.o - 


B,0 


= 0 , 


1 


1 


d% 


COSt/- 


d‘^, 


B.2 


d(j> 


- -h '^n _2 — 


<'3,2 


d% 


drj 


1 — eg 


cost 


B,k 


d(j) 


+ ^B,k - % 


B,k 


1 — 677 
1 

1 — 677 
1 

1 — 677 


■ cost 


3,0 


d9 


coscj)- 


d‘^, 


3,1 




3,0 


d9 


d% 


■ cost 


3,fc-l 


dt 

d% 


B,k-2 


d9 


dt 
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where 

(3.46) = ^^i^{t,r],e,(j))d(t). 

It is important to note the solution j. depends on e and this is the reason why we add the superscript e 
to C/|, and 

3.5. Initial-Boundary Layer Expansion. Above construction of initial layer and boundary layer yields 
an interesting fact that at the corner point (t,x) = (0,a;o) for Xq £ dfl, the initial layer starting from 
this point has a contribution on the boundary data, and the boundary layer starting from this point has a 
contribution on the initial data. Therefore, we have to find some additional functions to compensate these 
noises. The classical theory of asymptotic analysis requires the so-called initial-boundary layer, where both 
the temporal scaling and spacial scaling should be used simultaneously. Fortunately, based on our analysis, 
the improved compatibility condition (1X71) implies the value at this corner point is a constant for w ■ n < 0. 
Then These contribution must be zero at the zeroth order, i.e. 

(3.47) '^/o(r,fo,w) = 0, 

(3.48) = 0, 

Therefore, the zeroth order initial-boundary layer is absent. 


3.6. Construction of Asymptotic Expansion. The bridge between the interior solution, the initial layer, 
and the boundary layer is the initial and boundary condition of (HU. To avoid the introduction of higher 
order initial-boundary layer, we only require the zeroth order expansion of initial and boundary data be 
satisfied, i.e. we have 


(3.49) U^{0,x,w) + %^q{0,x,w) + ^gQ{0,x,w) 

(3.50) UQ{t,xo,w) + %^o{t,xo,w) + ‘WBfi{t,xo,w) 

The construction of C/|. and are as follows: 


Assume the cut-off function ip and ipo are defined as 


(3.51) 

(3.52) 

and define the force as 

(3.53) 


tPifi) = 


1 0<fi<l/2, 

0 3/4 < fj. < oo. 



0< fi< 1/4, 
3/8 < fi < oo. 


F{e]r]) = 


1 — £77 ’ 


h(x, w), 
g{t,xo,w). 


Step 1: Construction of zeroth order terms. 

The zeroth order boundary layer solution is defined as 


(3.54) 


< sin - F{e; rj) cos -h - fS 

/o^(LO,0,/-) 
lim^_,.oo/o(i, ??, 0, </) 


V’o(eb) - /o(Loo,6»)^, 

0 , 

for sin^ > 0, 
foit,oo,0). 


Assuming g £ L°°, by Theorem 14.11 we can show there exists a unique solution fQ{t,g,9,(j)) £ L°°. Hence, 
Q is well-defined. 
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The zeroth order initial layer is defined as 


(3.55) 


drH 

X, w) 
T^{0,x, w) 
^ \imr^ooJ^oiT,x,w) 


T^{t,x,w) - J^^{oo,x) 

0 , 

e"^Jo"(0,f, w) + (1 - e"^)Jo"(0,f), 
h{x, w), 

J^^{oo,x). 


Assuming h G L°°. Then we can show there exists a unique solution JFq{t,x,w) G L°°. Hence, '^/g is 
well-defined. 

Then we can define the zeroth order interior solution as 

C^o = U^, 

= 0 , 

x) = J^q(oo,x) in fl, 

Co) = foit,oo,0) on dn, 

where {t, x, w) is the same point as (r, rj, 0, (f). Note that due to the improved compatibility condition (I3.7L 
we have '^J_q( 0, ry, 0, 0) = '^/g(r, foj w) = 0. 


(3.56) 


dtU^ - A, 

^o(0, 


Step 2: Construction of first order terms. 

Define the first order boundary layer solution as 




df^ df^ 

(3.57) sincj)-^ - F{€]r]) cos + fl - fl 

flit, 0, 0,(1)) 

fl{t,r],0,(j)) 


'0o(e?7) yfi(t,ri,0,(l)) - /i"(t,oo,6») 
tjj^er]) d‘^B,o 


= COS((> 


1 — eTy 80 
w ■ VxUQ{t,XQ,w) for sin(/)>0, 
flit, 00,0). 


Define the first order initial layer as 


(3.58) < 


= J-liT,X,w) — Tl{oo,x) 


%\iT, X, w) 

drH 

F{iT,x,w) = e~'^ Fli0,x,w) + I — w ■yx'^ifi\is,x,w)e^~'^ds, 


w ■ Vx'^ffi I die, 


limT 


J^i(0, X, w) 
, 7'i(r, X, w) 


w ■ S/xUq{G,x, w), 
Ff{oo, x). 


Define the first order interior solution as 


(3.59) 


r _ ui 

I dtUi-^xUi 

1 U{id,x) 

I uiit,x) 


Ul-w- VxUS, 

0 , 

Tlico,x) in n, 
flit, 00 , 0 ) on 9f2. 


Step 3: Construction of ‘^2 ^^id 17|. 

Define the second order boundary layer solution as 

(3.60) 


'^B,2it,'n,^,4>) 

< sin - Fie; 7y) cos (j)-^ + /| - /| 

flit, 0, 0,(1)) 

\imr,^^flit,r],0,(j)) 


’fpoi^v) (^f2it,V,S,(l)) - flit,00,0)'^, 

ifierf) I ^^B,o 

COS(^- -- 

1 — ery 80 8t 

w ■'^xUiit,Xo,w) + 8tUlit,xo,w) for sintyi > 0, 
flit, 00,0). 
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Define the second order initial layer as 

= F^{t,x,w) - F^{oo,x) 

w ■ Six'll I )dw, 


(3.61) < 


drF^ = - 


J'Kr, X, w) = e J'KO, X, wJ) + j [ T 2 — w ■ x^^i.i ] {s, x,w)e^ "^ds, 


T 2 { 0 ,x,w) = w ■VxUl{0,x,w) + dtUQ{0,x,w), 

limT-_>oo-^2(D^, w) = 7)1(00, f). 

Define the first order interior solution as 

C/| = Ul-w-VxU^-dtUS, 

dtUi-AxU^ = 0, 

1 / 2 ( 0 , x) = 7 - 2 ( 00 , x) in fl, 

1 ) 2 ( 1 , x) = /|(t,oo,0) on dfl. 


(3.62) 


Step 4: Generalization to arbitrary k. 

Similar to above procedure, we can define the order boundary layer solution as 
(3.63) 


dft dft - i/'(er?) 9 '^r A- 1 t 2 

fl(t,0,9,(j)) = w-WxUl_^(t,xo,w) + dtU^_ 2 (t,xo,w) for sin()) > 0, 

lim^^oo/|(i, ly, 0) = fl(t,oo,e). 

Define the order initial layer as 


(3.64) < 


‘^i(k(T,x,w) = FI(t,x,w) - Tl(oo,x) 

w ■yx'^ik-A<)w, 


drn = - 


lim 


T^(t,x,w) = e '^T^(0,x,w) + J \^l — w-Vx‘77i^k_i](s,x,w)e’^ ’’ds 

jF^(0,x,w) = w-yxU^_i(0,x,w) + dtU^_ 2 ( 0 ,x,w), 

t^oo)^1(t,x,w) = Fl(oo,x). 


Define the k*^ order interior solution as 
(3.65) 


Ul = U^,-w-VxU^k-i-dtUI-2^ 


dtU^-AxUI = 0 , 

Ul(0,x) = Fl(oo,x) in ft, 

= /fe(^,oo,6») on on. 

When g and h are sufficiently smooth, then all the functions defined above are well-posed. The key point 
here is in the boundary layer, the source term including j, is in L°° due to the substitution (13.391) . 

4. e-MiLNE Problem 

In this section, we study the e-Milne problem for f’^(g,9,(j)) in the domain (g,9,(j)) G [0,oo) x [— 7r,7r) x 

[-7r,7r) 


(4.1) 

where 

(4.2) 


sincj)^ + F(e;g)cos(j)^ + - r = S'^(g,9,<j)), 

f'^(0,9,4>) = H'^(9,4>) for sin(^ > 0, 

lim^_>oo/"(??, 6*, <?i’) = fL(^), 
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(4.3) 

r./ ^ efj{er]) 

F{e;r]) = 

1 — er] 

(4.4) 

/ 1 0 < /r< 1 / 2 , 
1 0 3/4 < yi < 00 

(4.5) 

\H^{9,cf)\<M, 

and 


(4.6) 

\S^iri,9,cf)\<Me-^^, 


for M > 0 and K > 0 uniform in e and 9. In this section, since the key variables here are i] and cj), we 
temporarily ignore the dependence on e and 9. We define the norms in the space ( 77 , (jj) G [0, 00 ) x [—tt, tt) as 
follows: 

( 1^00 ^TT \ 1/2 

Jo J ’ 

(4-8) II/IIl-l- = sup |/(?7,^)l- 

{r],(j))^[0,oo) X [—7r,7r) 

In [131 Section 4], the authors proved the following results: 

Theorem 4.1. There exists a unique solution f{ri,(j)) to the e-Milne problem lEZF satisfying 

(4.9) \\f-fooh.L.<c(^l + M+^y 

Theorem 4.2. There exists a unique solution /(ly, <)>) to the e-Milne problem lEZP satisfying 

(4.10) \\f-U\L^L^<c(^l + M+^y 

Theorem 4.3. For Ko > 0 sufficiently small, the solution f(j],4>) to the e-Milne problem EIP satisfies 

(4.11) ||e^'>’'(/-/oo)|L.^. <C (^1 + M+^^, 

Theorem 4.4. For Kq > 0 sufficiently small, the solution f{r], (f) to the e-Milne problem o satisfies 

( 4 - 12 ) <c(^l + M+^y 

Theorem 4.5. The solution f{rj, (f) to the e-Milne problem B with S = 0 satisfies the maximum principle, 

i.e. 


(4.13) 


min h{(j)) < f(r],(j)) < max h{(j)). 

sin^>0 sin^>0 


Remark 4.6. Note that when F — 0, Theorem B Theorem \4.2[ Theorem B Theorem B and Theorem 
\4-5\ still hold. Hence, we can deduce the well-posedness, decay and maximum principle of the classical Milne 
problem 


(4.14) 


] /(O,0) 

I lim^^oo/(i?,())) 


h{(j)) for sin </> > 0 , 

/oo. 
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5. Diffusive Limit 

In this section, we prove the first part of Theorem 1 1.2 1 

Theorem 5.1. Assume g{t,xo,w) G (^^([OjOo) x T”) and h{x,w) G x 5^). Then for the unsteady 

neutron transport equation the unique solution u’^{t,x,w) € L°°([0,oo) x x S^) satisfies 

( 5 . 1 ) \\u^-u^-no-n,oh^ = oii), 

where the interior solution Uq is defined in US. 56\) . the initial layer is defined in i3.55\) . and the boundary 
layer q is defined in {3.54^ . 

Proof. We divide the proof into several steps: 


Step 1: Remainder definitions. 

We may rewrite the asymptotic expansion as follows: 


(5.2) 




The remainder can be defined as 

N 


/c=0 


N 


/c=0 


/c=0 


(5.3) 

Rn = u^-Y^e^Ut-Y^e^% 


O 

II 

o 

II 

where 


(5.4) 

Qn 

(5.5) 

^I.N 

(5.6) 

^B,N 


N 




N 


Qn = 

N 

N 


k=0 


Noting the equation is equivalent to the equations (13.211) and p.40|) . we write £ to denote the neutron 
transport operator as follows: 

Cu = e^dtu + ew ■ VxU + u — u 


(5.7) 

(5.8) 


= drUew ■ VxUu — u 


>du 


, du 


= e—+ sm^—-- 

at orj 1 — er] 


du du^ 


Step 2: Estimates of CQn- 

The interior contribution can be estimated as 

(5.9) £Qo = s^dtQo + ew ■ VxQo + Qo — Qo 

= e^dtU^ + ew ■ VxU^o + {U^ - Uq) = ™ 

We have 


(5.10) 

(5.11) 

This implies 

(5.12) 


\e^dtUS\ < Ce^\dtU^\<Ce^, 

lew-VxU^l < Ce\VxU^\<Ce. 


|£Qo| <Ce. 


Similarly, for higher order term, we can estimate 

(5.13) jrQN = e^dtQN + ew-yxQN + QN-QN = e^+^dtUf^ + e^+^w ■ VxU^n- 
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We have 

(5.14) 

(5.15) 

This implies 

(5.16) 


^N+l 


dtUU < Ce^+^\dtUl,\<Ce 


^N +2 


e'^+ize-V.C/^l < Ce^+^\V^U%\<Ce^+\ 
\CQn\ < Ce^+\ 


Step 3: Estimates of CJSi^n- 

The initial layer contribution can be estimated as 

(5.17) + ^/,o “ ^/,o 

= clr'^/o + ew ■ + '^/,o “ ^^1,0 = eVx'^/o- 

Based on the smoothness of '^/gj we have 

(5.18) |£j2/,o| = |eV,'^/o| < Ce. 

Similarly, we have 

(5.19) = dr^I,N + -Vx^I,N + ^I,N — ^I,N = 

Therefore, we have 

(5.20) \C^i,n\ = \e^^^Vx%\N\ < Ce^+^. 


Step 4: Estimates of C^b,n- 

The boundary layer solution is = (/^ — f^{oo)) ■ tpo = 'fk'fpo where f^{ri, 9, cj)) solves the e-Milne problem 

and '^k = fk~ /fc(oo)- Notice = V'o, so the boundary layer contribution can be estimated as 

(5.21) 


28 ^ 3,0 , . 
e —--h sm 


d^B 


dt 


2 ^ 
dt 
,d% 


f d^B,0 dJ^Bfi 

a, ■ --d —+ — 


1 _,, + -sri+ 


''0-^ I - --COS( 

orj / 1 — er] 


d(j) 89 

^- u , . J.I I d% , ^^8'ipo\ V'oV’e , 8% 

e — + smal j - — ^ 


■ V'o'^ — 


■ ^ 0 % — i’o'^o 


= e ~:yr + w simp— - 


8 t 

. 8 % 


eip 


^ 8 % 


Fi 1 C0S(()-—-b 1^-1^ -b sin 

ar] \ — erj op J orj 1 — erj 

. 8 % 


. 8 ^ 0 , 


V'oe , 8 % 

coaa^ 


2^'^0 , . idipo V'oe .oro 

e^ — -bsin(/>—ro - -cos^—. 

8 t or] 1 — er] 89 


It is easy to see 
(5.22) 


28 % 


8 % 

^ 8 t 

< 

8 t 


< Ce\ 


Since po = 1 when 77 < l/(4e), the effective region of 9,;V'o is 77 > l/(Ae) which is further and further from 
the origin as e —>■ 0. By Theorem 14.21 the first term in (15.211) can be controlled as 

, 8 ^ 0 . 


(5.23) 


s\np^—% 
8 r] 


< Ce- 


< Ce. 


Eor the second term in (j5.21L we have 

V'oe 


(5.24) 

This implies 

(5.25) 


1 — 677 89 


< Ce 

8 % 


89 


< Ce. 


\C^bA < Ce. 
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Similarly, for higher order term, we can estimate 


(5.26) C^B,N = 


dt 


■ sin(j)- 


= e 


N+2 


dt 


i=0 


drj 1 — erj 

dt] ' " 1 — er] 


. d^B,N d^B,N 
cosQi —-h 


d^N i ■ j.^^0 

+ > e sin(j)^—Ti - - -COS0 

nr) 1 — fr) 


d(j) do 

f 

~d 0 


It is obvious that 
(5.27) 


,A+2 

< e^+2 

d^N 

dt 


dt 


< Ce^+^. 


Away from the origin, the first term in (15.261) can be controlled as 

k 


(5.28) 


i=0 ' 


< Ce-^<Ce'^+\ 


For the second term in (|5.26|) . we have 

i/ioe'=+i 


(5.29) 

This implies 

(5.30) 


1 — er] dO 


< Ce'^+i 

dn 


dO 


< Ce'^+K 


\C^b,n\ < Ce'^+i. 


+ ^B,N — <Sb, 


N 


Step 5: Synthesis. 

In summary, since Cu’^ = 0, collecting (15.3p . ()5.16p . (j5.20L and (I5.30L we can prove 

(5.31) \CRN\<Ce^+\ 

Consider the asymptotic expansion to = 4, then the remainder R 4 satisfies the equation 

(5.32) 

edtR^ -\- ew • xRa T Ra: — 7^4 = /1774, 

i?4(0,f,w) = Efc=l 

Ri(t,XQ,w) = Y^\^i(^‘^ik{t,xo,w) for w-n<0 and xq G dfl. 

Note that the initial data and boundary data are nonzero due the contribution of initial layer and boundary 
data at the point {t,x) = (0,xo)- By Theorem 12.91 we have 


(5.33) l|774|lL~([o,oo)xnx5i) — ^'(fl) ||£i?4||j;^2([Q g^jxnxSQ + 

4 


' ([ 0 ,Oo) X ) 


k^l 

4 




k=l 


L“(nxsi) 


L°°([ 0 .i]xr-) 


< C{n) ( ^(Ce®) + (Ce®) ]+Ce + Ce = C{n)e. 


Hence, we have 

(5.34) 

Since it is easy to see 

(5.35) 




fc=0 


fc=0 




= o(l)- 




fc=l 


fc=l 


fc=l 


L°°([0,oo)xf2XtS^) 


= 0(e), 


L°°(nxSQ 

















































26 


LEI WU 


our result naturally follows. 


□ 


6. Counterexample for Classical Approach 

In this section, we present the classical approach in [1] to construct asymptotic expansion, especially 
the boundary layer expansion, and give a counterexample to show this method is problematic in unsteady 
equation. 

6.1. Discussion on Expansions except Boundary Layer. Basically, the expansions for interior solution 
and initial layer are identical to our method, so omit the details and only present the notation. We define 
the interior expansion as follows: 


k=0 


(6.1) U{t,x,w) ^ y^ e’^Ukit,x,w), 

Uo{t,x,w) satisfies the equation 

( 6 . 2 ) 

Ui (t, X, w) satisfies 
(6.3) 


Uq — Uq, 
dtUo-A,Uo = 0 . 


Ui = Ui-w-VxUo, 

dtUi-AxUi = 0 , 


Uk — Uk W'^xUk—1 dtUk — 2-: 
dtUk - AxUk = 0. 


and Uk{t,x,w) for fc > 2 satisfies 

(6.4) 

With the substitution (I3.19L we define the initial layer expansion as follows: 

OO 

(6.5) ^j(r,x, w) ^ y^ e''^i,k{T,x,w), 
where satisfies 

(6.6) 


k=0 


dr^I,0 = 0, 

^ i ^ o { t , x , w ) = e“'^'^/,o(0,f, w) + (1 - e“'^)‘^/,o(0,x). 

and ^i^kiT,x,w) for fc > 1 satisfies 


dr^i.k = — 


L 


w ■ Vx'^i,k -1 did 


^i^k{j,x,w) = e '^'^/,fe(0,f,rZ;) + / (- w • ) (s, x, w)e'’ '^ds. 


(6.7) 


6.2. Boundary Layer Expansion. By the idea in [I], the boundary layer expansion can be defined by 
introducing substitutions p.32p . (I3.34p . and (I3.37p . Note that we terminate here and do not further use 
substitution (13.391) . Hence, we have the transformed equation for (11.11) as 


( 6 . 8 ) 


e2_+sin(d + C) 


du'' e e 1 


= 0 , 


= g{6,^) for sin(d + ^) > 0. 


We now define the Milne expansion of boundary layer as follows: 

00 

^ (t, T ], 6», 0) - ^ e '"‘^ k { t , ??, 9, 4>), 




(6.9) 
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where can be determined by comparing the order of e via plugging (I6.9p into the equation (16.81) . Thus, 
in a neighborhood of the boundary, we have 


( 6 . 10 ) 

( 6 . 11 ) 

( 6 . 12 ) 

(6.13) 

where 


sin(6( + ^) 


dij 


+ ~ ^0 — 0, 


sin(6( + - '^1 

or] 

sin(6* + 

Of] 

sin(0 + C)^+‘%-'% 

or] 


1 

cos{0 + 4) 


1 — €77 
1 

1 — €77 
1 

1 — €77 


cos(0 + 4) 
cos(0 + 4) 


89 ’ 

d'^i 8% 
8‘^k-i d^k-2 


89 


8t 


(6.14) '^k{t,r],9) =-^ J '^fe(t,77,6»,4)d4. 

6.3. Classical Approach to Construct Asymptotic Expansion. Similarly, we require the zeroth order 
expansion of initial and boundary data be satisfied, i.e. we have 

(6.15) UoiO,x,w) +‘Wi^oiO,x,w) +'^B,o{0,x,w) = h, 

(6.16) Uo{t,xo,w) +‘^ifi(t,xo,'w) = g. 

The construction of Uk^ ^i,k, and '^B,k by the idea in [T] can be summarized as follows: 

Assume the cut-off function rlr and rjjo are defined as (13.5111 and (13.521) . 

Step 1: Construction of zeroth order terms. 

The zeroth order boundary layer solution is defined as 


(6.17) 


%(t,i],9,^) 

< sin(6»-h4)^-H/o-/o 

foit,0,9,0 

fo{t,r],9,0 


The zeroth order initial layer is defined as 


(6.18) 


‘Wifiir, X, w) 
8rTo 
Fair, X, w) 

Fo[0,x, w) 
lim^_>oo.?b(T, X, w) 


= 11:9,0 - fo{t,00,9)^, 

= 0 , 

= g{t,9,0 for sin(0-I-4) > 0, 

= fo{t,oo,9). 

Fo{t, X, w) - Fo{oo, x) 

0 , 

e”"” J'o(0, X, w) -t (1 - e""^)^^(0, 0: 

h{x, w), 

Fo{oo,x)- 


Then we can define the zeroth order interior solution as 


Uo 

8tUo — AxUo 

Uo{0:0 

Uoit,xo) 

where {t,x,w) is the same point as (r, 77 , 0 , 4 ). 



Uo: 

0 , 

Fo(oo,x) in n, 
fo{t,oo,9) on 8n, 


Step 2: Construction of first order terms. 
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Define the first order boundary layer solution as 


( 6 . 20 ) 


sm(0 + O—+/ 1-/1 - cos(0 + Oy3^^> 

= w ■VxUo{t,xo,w) for sin(6* + $) > 0, 

lim^^oo= fi{t,oo,e). 


Define the first order initial layer as 


'^i,i{j,x,w) = Ti{t,x,w) — Ti{<x>,x) 


drTi = - J ■Voo'^i^ojdw, 

J^i{t,x,w) = e~'^J^i{0,x,w)+ J — w ■Vx‘^i,o^{s,x,w)e‘^~'^ds, 

Ti{0,x,w) = w ■VxUq{0,x,w), 

liiui-^oo .Fi (r, X, w) = J'i(oo,x). 


Define the first order interior solution as 


Ui = Ui—W'VxUo, 

dtUi-AxUi = 0, 

17i(0,x) = J^i(oo,x) in fl, 

Ui{t,x) = fi{t,oo,0) on on. 


Step 3: Construction of second order terms. 

Define the second order boundary layer solution as 


= V’o(e»7) /2(t,?7,6»,^) -/2 (Coo,6») 


,,,3. .L-L = 

f2{t,O,0,^) = w-VxUi{t,xo,w)+dtUoit,xo,w) for sin(6» + ^)>0, 
lim^^oo= f 2 {t,oo, 0 ). 

Define the second order initial layer as 


^I, 2 {t,X,w) = T 2 {t,X,w) — 7 ^ 2 ( 00 , x) 


drJ^2 = -J \^W ■Vx‘^I,ljdw, 

J^2{t,x,w) = e~'^ 7^2(0, x,w) + [ (^2 — w ■'Vx‘^i,i^{s,x,w)e^~'^ds, 


7 ^ 2 ( 0 , X,w) = w • Va;17i(0,f,'u;) + 9tCfo(0,x, iZ;), 
limT-^oo-^2(D^, h;) = J'2(oo,x). 


Define the first order interior solution as 


(6.25) 


U 2 = U2-w-VxUi-dtUo, 

dtU2-AxU2 = 0, 

1 / 2 ( 0 , x) = 7 ^ 2 ( 00 , x) in n, 

U 2 (t,x) = f 2 (t,oo, 0 ) on i9D. 


Step 4: Generalization to arbitrary k. 
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Similar to above procedure, we can define the order boundary layer solution as 
(6.26) 

‘^k{t,ri,9,0 = ipo{ev)(9,0 - fk{t, 00 , 9 ) 


wfl , , r f /a , d^k-i d‘^k-2 

s,„(9 + 0 _ + A-/, = „s(9 + {)—- 

fk{t, 0,9,0 = W ■^xUk-i{t,xo,w) + dtUk-2{t,xo,w) for sin(6i + ^) > 0, 
fkit,V,9,0 = fk{t,oo,9). 

Define the order initial layer as 

‘^i^k{T,x,w) = Tk{T,x,w) - Fk{oo,x) 


(6.27) < 


drTk = - 


/51 


w ■ Vx'^Lk-i dw 


Tk{T,x,w) = e ^Tk{0,x,w) + J yJ^k - w ■Wx‘^i,k-i]{s,x,w)e'' '^ds 
^■^(O,^, u;) = w ■VxUl_j^{0,x,w) + dtUl_ 2 iO,x,w), 

\imr^ooJ^kiT,x,w) = Tk{oo,x)- 

Define the order interior solution as 


(6.28) 


Uk — Uk rC-Va;t//e —1 dtUk—2, 

dtUk — AxUk = 0, 

Uki0,x) = Fk{oo,x) in D, 

Uk = fk{t,oo,9) on dVt. 


By the idea in [T], we should be able to prove the following result: 

Theorem 6.1. Assume g(t,xo,w) and h(x,w) are sufficiently smooth. Then for the unsteady neutron 
transport equation a, the unique solution u^(t,x,w) € L°“([0, 00 ) x D x 5^) satisfies 


(6.29) 


— Uq — = 0{e). 


Similar to the analysis in na Section 2.2], considering a crucial observation that based on Remark |4.61 we 
know that the existence of solution fi requires 

(6.30) ^ (^fo{t, 77, 9,0 - fo{t, 00 , 9)^ e D°°([0, 00)^ x [-tt, tt) x [-tt, tt)). 

This in turn requires 

(6.31) ^ G L°°([0,oo)^ X [-7r,7r) x [-7r,7r)). 

On the other hand, as shown by the Appendix of [T3], we can show for specific g, it holds that d^ifo ^ 
L“([0, 00 )^ X [—7r,7r) x [— tt, tt)). Due to intrinsic singularity for (16.171) . this construction breaks down. 

6.4. Counterexample to Classical Approach. 

Theorem 6.2. If g{t,9,(j)) = Pe~*cos(j) and h{x,w) = 0, then there exists a C > 0 such that 

(6.32) ||u^-t/o-'^/,o-'^B,o|lioo >O>0 

when e is sufficiently small, where the interior solution Uq is defined in h6.19\) . the initial layer is defined 
in 116.18\} . and the boundary layer g is defined in 1)6.17^ . 

Proof. We divide the proof into several steps: 















30 


LEI WU 


Step 1: Basic settings. 

By (16.171) . the solution /o satisfies the Milne problem 


(6.33) 


sin(0 + ^)-^ +/o -/o 
lim^_>oo fo{t,r],0,^) 


0 , 

for sin(0 + ^)>O, 
fo{t,oo,0). 


For convenience of comparison, we make the substitution (p = 0 + ^ to obtain 


(6.34) 


sm + Jo- Jo 
or] 

fo{t,O,0,cj)) 

fo{t,r],0,(l)) 


0 , 

g{t,0,cj)) for sin(/)>0, 
/o(t,oo,6»). 


Assume the theorem is incorrect, i.e. 


(6.35) 


lim ||(17o + ^ 1,0 + ^B.,o) - {U^ + ^/,o + ^b,o)\\l^ = «. 


We can easily show the zeroth order initial layer '^b,o = '^J.o = 0 h{x,w) = 0. Since the boundary 

g{t, 9, (p) = t^e“‘ cos p independent of 0, by (16.1711 and (I3.54L it is obvious the limit of zeroth order boundary 
layer /o(t, oo,d) and fQ{t,oo,9) satisfy /o(t, oo,d) = C'i(<) and fQ{t,oo,9) = C 2 {t) for some constant C'i(t) 
and C 2 {t) independent of 0. By (16.1911 . (13.5611 and solution continuity of heat equation, we can derive the 
interior solutions are smooth and are close to constants Uq = C'i(t) and t/g = C 2 {t) in a neighborhood 0(e) 
of the boundary with difference 0(e). Hence, we may further derive in this neighborhood, 


(6.36) 


lim ||(/o(oo) A'^) - ifoioo) + = 0. 


For 0 < g < l/(2e), we have ipo = 1, which means fo = ■% + /o(oo) and /q = + /o(oo) in this 

neighborhood of the boundary. Define u = /o + 2, 17 = /g + 2 and G = g + 2 = t^e“‘ cos p + 2, then u{g, p) 
satisfies the equation 


(6.37) 


. , du 

smp——\- u — u 
ag 

u{Q,p) 

lim^_>oo u{g,p) 


0 , 

G{pi) for sin(^>0, 

2 + /o(oo). 


and U{g,p) satisfies the equation 

,dU 


,dU 


(6.38) 

Based on (I6.36|) . we have 


sin p^ + F(e; g) cos p-— + U -U 
og op 

Ui0,P) 


lim^_),oo U{g, p) 


0 , 

G{p) for sin^ > 0, 

2 + /o(oo). 


(6.39) 


lini ||t7(ry,<))) - u{g,p)\\^^ = 0. 


Then it naturally implies 

(6.40) lim||c7(?7)-u(ry)|| =0. 


Step 2: Continuity of u and U at g = 0. 

For the problem (I6.37|) . we have for any rg > 0 

(6.41) |u(?7)-u(0)| < [ \u{g,p)-u{0,p)\dp+ [ \u{g,p) - u{0,p)\dp\ 

\Jsin<p<ro Jsin<p>rQ / 

Since we have shown u G L°°([0, oo) x [—tt, tt)), then for any d > 0, we can take rg sufficiently small such 
that 

1 f c s 

(6.42) — \u{g,p)-u{0,p)\dp < —arcsinrg<-. 

Jsin<p<ro ^ 
















GEOMETRIC CORRECTION FOR DIFFUSIVE EXPANSION OF NEUTRON TRANSPORT EQUATION 


31 


For fixed tq satisfying above requirement, we estimate the integral on sin^ > tq. By Ukai’s trace theorem, 
is well-defined in the domain sin^ > rp and is continuous. Also, by consider the relation 


(6.43) 


m(0) - U(0, 0) 

^(0,0) = -T-T-> 

orj sm (j) 


we can obtain in this domain drjU is bounded, which further implies u(? 7 , 0 ) is uniformly continuous at 17 = 0 . 
Then there exists Jp > 0 sufficiently small, such that for any 0 < t] < Sq, we have 

(6.44) ^ / |m( 77,0)-M(O,0)|d0 < -^ [ ^d0<^- 

Jsin<p>ro Jsin<p>rQ ^ ^ 

In summary, we have shown for any (5 > 0, there exists Jp > 0 such that for any 0 < ry < Jqj 


(6.45) 


Hv) 


«(0)|<^ + 


6 

2 


<5. 


Hence, u{ri) is continuous at 77 = 0. By a similar argument along the characteristics, we can show U{r], 0) is 
also continuous at 77 = 0 . 

In the following, by the continuity, we assume for arbitrary d > 0, there exists a <5p > 0 such that for any 
0 < 77 < 5o, we have 


(6.46) 177 ( 77 ) — u(0)| < S, 

(6.47) |C7(77)-t7(0)| < 6 . 


Step 3: Milne formulation. 

We consider the solution at a specific point ( 77 , 0) = (ne, e) for some fixed tt- > 0. The solution along the 
characteristics can be rewritten as follows: 

pne 1 

(6.48) u{ne,e) = — u{K)dK, 

Jo sm e 

ne 1 ne 1 1 

(6.49) U{n€,e) = G(€o)e~ -__i7(^)ci«:, 

Jo sm 0 (k) 

where we have the conserved energy along the characteristics 

(6.50) 75 ( 77 ,0) = cos 0e“^^’^\ 

in which (0, ep) and (C, 0(C)) are in the same characteristics of (ne, e). 


Step 4: Estimates of (16.481) . 

We turn to the Milne problem for u. We have the natural estimate 


(6.51) 


-{nt-K) 


1 


■d«: = 


sine 


= e 


= e 


nrif. 1 

/ e-r('^'^-«)-dK-ho(e) 

Jo e 

/ e“-dK-|-o(e) 

Jo £ 

/ e^dC -I- o(e) 

Jo 


= (I-e-") + o(e). 

Then for 0 < e < dp, we have | 7 t( 0 ) — u{k)\ < S, which implies 


(6.52) 


^ne 1 

e-snr7("'^-«)_ 


sine 


■u{K)dK = 


fUt 1 


sine 


■77(0)d«: -I- 0[5) 


= (1 - e-’^)7i(0) + o(e) + 0(d). 


For the boundary data term, it is easy to see 

(6.53) G(e)e-^”'= = e-”G(e)-k o(e) 

In summary, we have 

(6.54) n(ne, e) = (I — e“")M(0) + e“"G(e) -I- o(e) -I- 0{5). 
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Step 5: Estimates of (16.491) . 

We consider the e-Milne problem for U. For e << 1 sufficiently small, = 1- Then we may estimate 

cos(/)(C)e“^^^^ = cosee“^^"''^\ 


(6.55) 

which implies 

(6.56) 
and hence 


1 - ne^ 

cos0(Q = — -^cose. 


(6.57) 


• R - T1J7\ /e(ne-C)(2-eC-ne2) 

sin (()(0 = V1 - cos^ (/'(C) = W-- 

For C G [0, e] and ne sufficiently small, by Taylor’s expansion, we have 

(6.58) 1-eC = l + o(e), 

(6.59) 2-eC-ne^ = 2 + o{e), 

(6.60) sin^e = e'^ + oie^), 

(6.61) cos^ e = 1 — e^-I-o(e^). 

Hence, we have 


cos"^ e -I- sin e. 


(6.62) sin(/i(C) = ■\/e(e -|- 2 ne — 2Q + o(e^). 

Since y'e(e -I- 2ne — 2C) = 0(e), we can further estimate 

1 1 


(6.63) 


(6.64) 


sm 


1 


4>(C) \Je(e -I- 2ne - 2C) 


sin (/'(C) 


dC = 


e -I- 2ne — 2C 


+ 0 ( 1 ) 


+ o(e) = 1 - 


e -|- 2 ne — 2 k 


Then we can easily derive the integral estimate 


(6.65) 


Q Ik. sin d>(C) - 


1 


sin (/'(«') 


dre = e^ 


_ /.+2ne-2j 
^ V e 


ie^ 

Al+2n)t 

/ eT - 

2 . 

A 

ie^ 

^l+2n -| 

2 . 

A vp 


\/ e(e + 2ne — 2 k) 
1 


+ o(e). 


d/c -|- o(e) 


€<J 


da -\- o(e) 


/ •\/l+2n 

e“‘dt -I- o(e) 


= e 


= (l_ei-vT+^) + o(e). 

Then for 0 < e < (5o, we have \U{0) — U{k) \ < S, which implies 


( 6 . 66 ) 


n fn sliViTcfT _ 


1 




_ p 


ydC. 


1 


sinC'(K) 


sin(/)(K) 


{7(0)dK -I- 0(6) 


= (l_ei-+i+^)t/(0) + o(e)+O((5). 
For the boundary data term, since G((/') is G^, a similar argument shows 

(6.67) G(eo)e~'^‘' = e^“+^'''^"G(-\/l -|- 2ne) -I- o(e). 

Therefore, we have 

(6.68) Uine, e) = (1 - e^-^^)U(0) + e^-^^G(VT+^e) + o(e) + 0(5). 
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Step 6: Contradiction. 

In summary, we have the estimate 

(6.69) u{ne,e) = (1 — e“")u(0) + e“"G'(e) + o(e) + 0((5), 

(6.70) U{ne,e) = (1 - + 2ne) + o{e) + 0{6). 

The boundary data is G = t^e~* coscj) + 2. Fix t = 1. Then by the maximum principle in Theorem 14.31 we 
can achieve 1 < u(0, </>) < 3 and 1 < U{0,cj)) < 3. Since 


(6.71) 


w(0) = 


1 

27T 

_ 1 
27T 

= 2+^ 


w(0,= 


27r 


sin ^>0 


?x(0, (l))d(j) - 


27r 


sin <p<.0 


u(0, (^)d4' 


sin ^>0 
1 

27r 


(2 + e ^ cos 4>)d(j) 




sin (p>0 

1 


e ^ cos (/>d^ + 7 ^ 
2?! 


/ u{0,(j))d(j) 

sin <j5><0 

f u{0,(j))d(j), 

sin <^<0 


1 


2 — ^e ^ < u(0) < 2 + -e 


2-ie-i<i7(0)<2+ie-i. 


we naturally obtain 

(6.72) 

Similarly, we can obtain 

(6.73) 

Furthermore, for e sufficiently small, we have 

(6.74) G(V1 + 2ne) = 2 + e-i + o(e), 

(6.75) G(e) = 2 + e"i + o(e). 

Hence, we can obtain 

(6.76) u(ne,e) = ■u(O) + e“"(—tt(0) + 2 + e“^) + o(e) + 0(5), 

(6.77) 


U{ne,e) = U{0) + + 2 + e-^) + oie) + 0{S). 


Then we can see lime_>o ||G(0) — 'S(0)||^„„ = 0 naturally leads to lime_j.o ||(—?2(0) + 2 + e“^) — (—17(0) + 2 + e“^) 
0. Also, we have —u(0) + 2 + e“^ = 0(1) and —0(0) + 2 + e“^ = 0(1). Due to the smallness of e and 5, and 
also e“" ^ gi-Vi+2ra^ obtain 

(6.78) |0(ne, e) — u(ne, e)| = 0(1). 

However, above result contradicts our assumption that lime_>o \\U{r], </>) — 11 ( 77 , (I))\\l'=° ~ ^ '/')■ This 

completes the proof. 


□ 
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